Slip and flow of hard-sphere cofloidal glasses 
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We study the flow of concentrated hard-sphere coUoidal suspensions along smooth, non-stick 
walls using cone-plate rheometry and simultaneous confocal microscopy. In the glass regime, the 
global flow shows a transition from Herschel-Bulkley behavior at large shear rate to a characteristic 
Bingham slip response at small rates, absent for ergodic colloidal fluids. Imaging reveals both the 
'solid' microstructure during full slip and the local nature of the 'slip to shear' transition. Both the 
local and global flow are described by a phenomenological model, and the associated Bingham slip 
parameters exhibit characteristic scaling with size and concentration of the hard spheres. 

PACS numbers: 82.70.y, 83.50.-v, 83.60.-a, 83.85.Ei 



Wall slip in fluid flow has received considerable atten- 
tion for many years [E i, i i, i, i, S i, & [13, El 113, [11 • 
Even for simple fluids it has been realized that the no- 
slip boundary condition may be inappropriate on length 
scales relevant for nanoporous media or nanofluidics [l|. 
More widespread is the presence of slip in complex fluids 
like suspensions and emulsions |3. Isl. B. la. m. YA Isl I Ol. I lOl I 111. 
1% Il3j . Despite the more accessible length scales in these 



systems, it has remained challenging to gain microscopic 
insight into the nature of slip and understand its depen- 
dence on material composition, wall properties and flow 
rate. Recently, considerable progress has been made for 
soft particle pastes [7J, but for hard sphere (HS) suspen- 
sions the situation remains unclear. Slip was observed 
both in solid- llillJ] and liquid-like 13] particulate sus- 



pensions, but its possible relation to Brownian motion 
and the glass transition was mostly ignored. 

Experimentally, a proper interpretation of slip requires 
not only precise rheological data but also detailed spatial 
characterization of the flow. Various imagirig methods 
have been employed the last decades [5|, |7|, [iSJ, |lj, |l5|, |l6| ; 
however they either lack the resolution to study flow near 
the wall on the single particle level [5|, [tI, [ij, [ij , or they 
do not provide simultaneous rheological and microscopic 
information on the same sample |15l . llq . 

Here, we use confocal imaging in a cone-plate rheome- 
ter to address slip and yielding of dense HS colloids on 
both the microscopic and macroscopic scales. Slip ap- 
pears as a yield stress emerges on entering the glass 
regime, with plug-flow persisting down to the colloid- 
wall interface. We flnd a 'Bingham' slip response, unlike 
emulsions JJ. The physics of slip is thus not universal 
between different classes of soft glassy materials. 

We used polymethylmethacrylate colloids of various 
radii (a = 138 nm, 150 nm, 302 nm, polydispersity 
~ 15%, all from light scattering), stabilized with poly- 
12-hydroxystearic acid, suspended in a mixture of decalin 
and tetralin (viscosity rjs = 2.3 mPas) for refractive index 



(RI) matching and seeded with ^ 0.5% of fluorescent col- 
loids of the same kind (a — 652 nm) for confocal imaging. 
In these solvents, the colloid interaction is very nearly 
HS-like [17] . Batches of different volume fractions (j) were 
prepared by diluting samples centrifuged to a random 
close packed sediment, taken to be at 0rcp — 0.67 J18J . 
When comparing data from different colloid sizes, how- 
ever, we report results in the reduced variable 1 — 0/(/),cp, 
which is independent of the numerical value of (/),cp- 

Measurements were performed in a controlled stress 
rheometer (AR2000, TA Instr.) in cone-plate geometry 
(radius Vc ~ 20 mm, cone angle 9=1°) with a mod- 
ifled base on which a glass slide (radius 25 mm, thick- 
ness ~ 180 /im 19], local roughness < 1 nm from AFM) 
is mounted. A solvent trap minimizes evaporation. By 
coupling in a piezo-mounted objective (> 60x, oil im- 
mersion) and optics via an adjustable arm connected to 
a confocal scanner (VT-Eye, Visitech Int.) we measure 
the velocity profile v{z,r) (with z the velocity-gradient 
direction) from movies taken at 5z = 2-5 /im intervals 
at a frame rate < 90 Hz at various distances r from the 
cone center, Fig.[lja). To prevent slip, both the glass and 
the cone can be made rough on the colloid scale by spin- 
coating a, (p r^ 0.3 suspension and sintering the resulting 
disordered colloidal monolayer. Experiments with both 
smooth and coated glass plates were performed; the cone 
is always coated to ensure stick boundary conditions at 
the top. All data were collected with controlled applied 
shear rate 7a (going from high to low 7^) but stress con- 
trolled measurements gave the same results. 

We first discuss how the rheology depends on and the 
wall conditions. Fig. (^b) presents the measured stress 
Gm versus 7a for coated and uncoated glass slides at vari- 
ous (/) for a = 138 nm. For a concentrated fluid ((/> = 0.52) 
below the colloidal glass transition {4>g ~ 0.57 from mean 
squared displacement measurements) we find linear be- 
havior at the smallest 7a and shear thinning at higher 
7a, independent of boundary condition. 




FIG. 1: (a) Cone-plate rheometer with transparent base 
and optics connected via an adjustable arm to the confocal 
scanner, (b) Measured flow curves ami'ja) for colloids with 
a = 138 nm, at various (j) for coated (Bj^jA) and un-coated 
plates (n,o,A). Dashed line: viscous flow. Full line: am from 
Eq. (|4]) for (j) — 0.59, using Eq. ((5)1 for r > ^^(71,) and param- 
eters as = 1.8 Pa, /3 = 8.2 • lO"* Pa s m"\ ay = 5.5 Pa and 
a = 6.1 Pa s^". Dotted line: critical applied rate ^a,c (see 
text), (c) Reduced stress {am — as)/rjcs versus 7a for (f) = 0.58 
(V), 0.59 (o), and 0.64 (D). Line: {a,^ - as)/VcS = >. 



> (f)g the coated and uncoated results 

With coated glass and cone, we ob- 

13l), Hcrschel-Bulkley (HB) behavior: 

with cr„ the yield stress. For > 61%, 

< 1.2, to be 



However, for 
differ markedly, 
serve (as before 

Cm == o'l/ + Q!7a 1 with ay the yield stress, 
shear localization occurs for 1 < cr,n/cy 



discussed elsewhere [21[. On the other hand, for the 
smooth surface, cTm^ja) shows two regimes. For stresses 
somewhat below ay, we find apparent flow described by 
Cm = cTs + rjeS'Ja, Fig. [BJc), with a threshold as and 
an effective viscosity rjes (regime I). As explained below, 
this corresponds to full wall-slip and solid-body rotation 
of the sample over the entire geometry. This 'Bingham' 
slip differs from non-Brownian suspensions |12l | , where no 
threshold as is seen. It is also distinct from the behav- 
ior a„i — cTg ex y/j^ in soft particle pastes [7|. In regime 
II, arnija) dcviatcs from the Bingham form, the sample 
starts to yield and a^n approaches the HB curve. 

To check if the slip in regime I reflects full sliding from 
the first colloid layer or, instead, local yielding of the 
structure near the wall, we used the fluorescent batch 
(a = 652 nm, (j) = 0.63) and performed 3D imaging and 
tracking [24I of the fiow directly above the glass. The 
density profile. Fig. [2lja), shows that the surface induces 
layering. Nevertheless, the corresponding microscopic ve- 
locity profile, Fig. [2l^a) , reveals full sliding with respect 
to the glass, with v ~ Vcone down to the first layer and 
no yielding. Further analysis (data not shown) confirmed 
that colloids remained caged. The data are unchanged 
for an RI and density matching solvent, i.e. sedimenta- 
tion is not affecting these larger particles. Moreover, as 
the large and small particles display the same Bingham 
response in Regime I (data not shown), we expect that 
the behavior in Fig. [2] is particle size independent. 



We now turn to the velocity profiles associated with 
the rheology in Fig. [Tfb). Fig. [2jb) presents v{z)/vcoaa 
at a distance r = 2.5 mm from the center of the cone for 
(f) = 0.59 at various 7a. We observe linear profiles for all 
7a but with a finite slip velocity, defined as the z = 
intercept of v{z), and a slope corresponding to a bulk 
shear rate 7 < 7a. On reducing 7a, the profiles exhibit 
increasing slip and eventually show a transition to plug 
flow, 7 = 0, v{z) = Wconc, as in Fig. Wi,&)- For 4> > 61%, 
around the transition, we observe shear localized near the 
cone, while away from the cone we have plug flow with 
w = Us > 0.7uconc down to the plate. Somewhat below 
ay, we again find Vg = v{z) = Ucono, sec Fig. [2ja). 

As seen for the profile at 7a = 1 s^^ in Fig. [2l^b), plug 
flow can occur locally even when the global flow curve is 
close to that for rough walls. To bring out the behavior 
as function of r when entering regime II, we show in 
Fig. El^a) the local shear rate obtained from imaging at 
various distances r at fixed 7a — 1.1 s~^. Solid body 
rotation (7 = 0) is present at small r, while 7 > for 
larger r. This evidences that, at fixed 7a, the stress is 
r-dependent [3| , even where the sample has yielded. 

The local and global rheology described above can 
be rationalized using a simple model similar to [3|, Il2l |. 
First, as suggested by the phenomenology of the data 
in Figs. [ljb),(c), we relate the local stress a to the slip 
velocity Vs of the colloids along the wall by: 



a = as + Pvs, 



(1) 



with a hydrodynamic term (3vs and a threshold stress as 
For the bulk flow, we use the HB form with n = 0.5: 



a = ay + a-^ 



0.5 



(2) 
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FIG. 2: (a) Density profile P(z) (D) during slip, from 3D- 
imaging of the a — 652 nm system ai tf) — 0.63, r = 2.5 mm 
and 7a — 0.01 s~^. (•) Corresponding velocity profile, show- 
ing full plug flow. Inset: Slice of a 3D-image, showing the 
flrst colloid layer. Scale bar: 10 jj,m. (b) Velocity profiles for 
the (j) = 0.59 data in Fig. [Hb), in units of the cone velocity 
Vcone = jaOr, as functlon of reduced height z/zg^p = z/9r, for 
various applied rates 7a at r = 2.5 mm. The arrow marks the 
slip velocity for 7a = 30 s~^. Dashed line: behavior without 
slip V = 7aZ as observed for (j> < (j)g ov in newtonian liquids. 
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FIG. 3: (a) Local shear rate 7, obtained from the velocity 
profiles, versus position r for 7a = 1.1 s~^, a = 138 nm and 
(f} = 0.59, the line is calculated from Eqs. Q,®. (b) Reduced 
local shear rate F = 2{/36r/a)'^'y versus reduced applied rate 
Q, = 2{l36/a)^r{r — ry)'ya in regime II, for two colloid sizes a, 
various and distances r. Included are data at fixed ja but 
different r. Line: the prediction F = 1 + f2 — Vl + 2f2. 



The local bulk shear rate j{r) and Vs{r) are related by 



iir) = ja 



6r 



(3) 



We solve the flow by approximating, at each r, the cone- 
plate as parallel plates with a gap Or and balancing the 
stress in the bulk and at the wall. For CTs < o'(r) < ay, 
7(r) =0 and the stress is given by Eq. H]). Shear starts 
when the local stress cr(r) induced by slip exceeds ay, 
i.e., when Vs > Vs = {ay — as)/ /3, which in a cone-plate 
is equivalent to r > rj^ = {ay — as)/{(3jad)- Shp and 
shear are then both present, and balancing the stress in 
Eqs. (HI),© gives Vg. The measured stress am is thus: 



- Pvs{r)]2rdr, 



(4) 



with Vs{r) — 'jaOr for r < ry, while for r > ry we have: 



Vs{r) = iadry 



2P^0r 



4 . {r 



.(5) 



The transition from regime I to II occurs when ry 



i.e. for ja,c = "ln!!° and a^ = 



2(Ty +crs 



. Eq. ^ also gives 



(iBr^ """" - m — 3 

the relation ri^a = ^fiOrd'^- A calculation of the radial 
shear ^9^7 reveals that the parallel plate approximation 
is valid for all r except ry < r <ry{\ ^^ 26). 

To test this model, we extracted the parameters (/3, as) 
and (a, ay) from the smooth and rough wall flow curves 
respectively, for various (j) and a. Figure |4] shows the (j) 
dependence of {(3, as), which we discuss later. We veri- 
fled the form ryeS = 2/30rc/3 using different cones. The 
predictions for the smooth- wall flow curves, using Eq. (|4]) 
and the extracted parameters, agree well with measure- 
ments; an example is shown in Fig. [T]J^b). Additionally, 
Eqs. Q,® reasonably describe the r-dependence of the 
local shear rate 7 in Fig. Wia); the predicted curve fol- 
lows directly from the rheological parameters, without 



fitting. We also test the dependence of the local shear 
rate 7 on 7a by plotting the reduced local shear rate 
r = 2(/36'r/a)^7 versus the reduced applied rate fi = 
2{l39 / a)'^r{r - ry)'ja, Fig-Eljb). With this normalization, 
our model predicts a master curve F = 1 -f Jl — ^/l + 2fl. 
The data indeed follow this behavior for various (j) and r, 
lending strong support for our model. For (j> > 0.61 dis- 
crepancies can occur for smaller $7, e.g. due to localized 
shear near the cone around yielding 2l|, |23| . 



Summarizing thus far, a HS colloidal suspension driven 
along a smooth wall shows slip for (p > (pg, where a finite 
yield stress appears. The global and local rheology at 
different concentrations and colloid sizes in a cone-plate 
geometry are described by a model with a Bingham re- 
lation between the local stress and slip velocity, which 
should be applicable to flow of HS's in other geometries. 
As evidenced in Fig. [21 the Bingham form is associated 
with plug flow down to the flrst colloid layer. 

The slip parameters (3 and as in Fig. [4] both increase 
strongly for cj) — > (/)rcp, but their physical origins are less 
clear. A possible source for as is van der Waals attraction 
2J| , but we minimized this by RI matching. Experiments 
with an RI mismatch 5n < 0.01, controlled by changing 
temperature {T), displayed no change in as either }25j . 
However, using decalin as solvent {Sn > 0.02), we ob- 
served no slip and a,n{ja) f^oUows the HB form at all ja 
for (j) > (j)g [2Q] . as with coated walls. Indeed, imaging di- 
lute suspensions with a = 652 nm at equilibrium showed 
that in decalin colloids were stuck to the glass, while with 
RI matching such sticking was absent. Therefore we must 
seek a different origin for as in our system. 

In an equilibrium HS fluid next to a wall, the contact 
value of the distribution function, n{0), is proportional 
to the pressure [2a]. For quiescent HS glasses next to 
walls, little information exists 27[. Nevertheless, con- 
sistent with Fig. [IJa), we expect colloid- wall contacts 
in our system. Therefore as may reflect a Coulomb- 
like friction associated with these contacts. The os- 
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FIG. 4: (a) Normalized slip coefficient Pa/rjs versus reduced 
volume fraction l — {(j)/(j}icp) for the different colloid sizes. Full 
line: Pa/ris — 0.9(1 — (j}/(j}^cp)~^ ■ (b) Normalized threshold 
stress asO? /ksT, for different colloid sizes (symbols as in (a)). 
Dotted line: aa oc 11, continuous line as oc H^'^^. 



motic pressure, Il{4>), may then play the role of 'normal 
force' in the Coulomb law, so we expect cts oc II. For 
HS glasses, II is uncertain, but a widely used form is 
n = 2.9no/(l - (t)/(l)rcp) [ii], with Ho = 3#B 7^/(4^3). 
Thus asa^/ksT should collapse data for different colloid 
sizes, as is indeed the case. Fig. |31Jb). But the data evi- 
dence as/Ho ~ (1 — 4>/4>rc-p)^^ with V ~ 2-3, apparently 
inconsistent with Coulomb friction. However, our col- 
loids are polydisperse, and the colloidal glass is under 
finite shear strain during slip, so that the above form of 
11(0) may be inappropriate. 

The parameter /3 reflects lubrication between the first 
layer of colloids and the wall. While these colloids fluc- 
tuate (with a fraction ~ n(0) in wall contact), we can 
deduce the thickness ^ of an effective lubrication layer 
via ^ — rjs/fi |12l |. Indeed, temperature-dependent ex- 
periments conflrmed the scaling j3{T) ex rjsiT) for fixed 
(j). As shown in Fig. \M.,a), scaling ^ by a collapses data 
(plotted as £,~^ oc 0) for different colloid sizes |29|, with ^ 
a small and decreasing fraction of a as </> — *■ (/)rcp • Empir- 
ically we find a/C ^ 0.9(1 - </)/0rcp)~^ ^ n/(3.2no). In 
principle, /? should follow from integrating the distribu- 
tion of colloid- wall gap sizes, n(z), with an appropriate 
form for the lubrication force [30|. However, this gives a 
logarithmic form j3a/ris oc — ln(l — (/i/(/),cp), in contrast to 
the data. Instead, we can identify ^ with the mean spac- 
ing between the first layer of colloids and the wall. The 
latter is obtained by assuming that the contact value of 
the density n(0) ex H and using for the density away from 
the wall, n{z)^ a form suggested in [26|. This equilibrium 
estimate gives a/^ ~ [</)(! — (j)/(j),-cp)]~^, very similar to 
our empirical form and to previous observations [13]. 

The above discussion of Cs and /? led to correct scaling 
of the data for different particle sizes. However, we did 
not obtain the correct functional dependence of cr^ on (j), 
and having to use a m,ean spacing to predict /3 appears 
rather ad hoc. Possibly other aspects of the physics must 
be taken into account. For example, we only observed 
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SO that the finite strain in the slipping 



sample due to 'cage elasticity l3l|] may play a crucial part. 

To conclude, we have demonstrated Bingham-like slip 
in dense HS suspensions above the glass transition, flow- 
ing along smooth, non-stick walls. The apparent simplic- 
ity of Eq. [TJ however, emerges from the complex physics 
of HS glasses next to walls, which is not yet well under- 
stood. Our findings contrast with observations in dense 
emulsions 7], showing that details matter in accounting 
for slip in soft glassy materials. Interestingly, emulsions 
may be HS-like at (/> < 0rcp [32], implying a possible 
crossover in slip response from that found in 7] to that 
reported here on reducing (j) in dense emulsions. 
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